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Abstract 

We study the supersymmetric circular Wilson loops of TV = 4 super Yang-Mills in large 
representations of the gauge group. In particular, we obtain the spectral curves of the 
matrix model which captures the expectation value of the loops. These spectral curves 
are then proven to be precisely the hyperelliptic surfaces that characterize the bubbling 
solutions dual to the Wilson loops, thus yielding an example of a geometry emerging 
from an eigenvalue distribution. We finally discuss the Wilson loop expectation value 
from the matrix model and from supergravity. 
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1 Introduction 

In gauge theory/gravity correspondences, a saddle point in the gauge theory path 
integral is expected to represent the space-time geometry in gravity. Since the saddle 
point is determined by the dynamics of the gauge theory, the space-time is said to be 
emergent. A notable example of such a phenomenon is the emergence of the sphere 
of the dual AdS^ x geometry from the matrix quantum mechanics governing the 
strong coupling dynamics of the constant modes of the scalars of A/" = 4 super Yang- 
Mills compactified on a S'^ [l]Jl] 

^See [2] for a review of subsequent developments and a list of relevant references. 
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When an operator is inserted in the gauge theory path integral, the saddle point, as 
well as the space-time represented by it, gets deformed. The new space-time develops 
bubbles of new cycles carrying fluxes. Such bubbling geometries were originally found 
for half-BPS local operators in A/" = 4 super Yang-Mills theory in the context of 
the AdS/CFT correspondence [3i|4j. They were later generalized to include Wilson 
loops [5||6|I7] and surface operators [8] of the Af = 4 theory, while bubbling in topological 
string theory was found and studied in p[TO|[TT]. 

The current work revisits the bubbling geometries for circular supersymmetric Wil- 
son loops in A/" = 4 super Yang-Mills. These geometries were constructed in a complete 
form in reference [7]. The ten-dimensional space-time is a warped product 

ds' = flds\,s, + fldsl, + fldsl, + dsl (1) 

of AdS2 X S"^ X and a half-plane S. The radii /i, /2, and all other supergravity 
fields are functions on E given in terms of two holomorphic functions, A and B. In 
fact, S is naturally identified with the lower half-plane in one sheet of a hyperelliptic 
surface, also denoted by S, and A and B are constructed geometrically. Thus the data 
(S, A, B) completely characterize the bubbling solution. 

In this paper, we demonstrate that the deformed saddle points in gauge theory 
represent the bubbling geometries by making use of a matrix model. It was conjectured 
in p^ [T3] that the Wilson loop expectation value is captured by the Gaussian matrix 
model with a loop operator insertion. The conjecture was recently proved in reference 
|14] . where it was also shown that the matrix is the constant mode of a scalar fieldEI 
We show that the saddle point configuration of the matrix eigenvalues back-reacts to 
the operator insertion and the hyperelliptic surface S arises as the spectral curve in a 
generalized sense that we explain in detaili We also find an interpretation of A and 
B in the matrix model. 

Concretely, the circular supersymmetric Wilson loop is defined as 

Wr = Tr^Pexp ^{iA + 9'4>'ds). (2) 

Here A is the gauge field and 0* are the six real scalars. The integral is along a circle 
in M^, 6^ is a constant unit vector in M^, and s is the parameter of the circle such 
that ||(ia;/(is|| = 1. The trace is taken in an irreducible representation R of U{N) or 
SU{N). Such R is specified by a Young tableau, which is also denoted by the same 

^ In [K] it was argued that the matrix model arises as a mirror of the topological A-model for the 
AdSs X superstring p^ . 

It was originally argued by Yamaguchi ^ that the eigenvalue distribution of the matrix model 
characterizes the bubbling geometry. 
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symbol R. The dual bubbling geometry has small curvature when R has long edges 
and it is characterized by a genus g hyperelliptic surface S, where g is the number of 
blocks in R (see fig. [1]). The Wilson loop expectation value is given by the matrix 
integral 

{WR)yM = [rfM]e-^^^^^Tr^e^. (3) 

The N X N matrix M is hermitian and Z is the partition function. For representations 
R that give rise to smooth bubbling geometries, we solve the matrix model in the limit 
where N is infinite and the 't Hooft coupling A = gyM^ is large. As it turns out, 
A and B are simply related to the resolvent uj{z) and the spectral parameter z of the 
matrix model: 

A oc uj{z) — 2z, B (X z + const. (4) 

We also show that the resolvent is given by the indefinite integral of a meromorphic 
1-form a on the same hyperelliptic surface S. The surface S is given by the equation 

= H2g+2{z) , (5) 

and the 1-form a by 



a = du;= 2-2 \ dz . (6) 




The polynomials H{z) and a{z) have degrees 2g + 2 and g + 1 respectively. We find 
from the matrix model analysis a set of constraints that determine the coefficients of 
a{z) and H{z) uniquely. These constraints are identical to the ones that arise in the 
bubbling geometry. The surface S is the spectral curve of the matrix model in the 
sense that the eigenvalue distribution is determined by S and a. 

Given our large solution of the matrix model, the Wilson loop expectation value 
can be easily computed. A natural question is whether it can also be reproduced in 
supergravity, by evaluating the on-shell action in the bubbling geometry background. 
We include in this paper some relevant calculations that will be useful for this purpose. 
In particular, we show that the on-shell supergravity Lagrangian is always a total 
derivative. This would imply that the on-shell action splits into two contributions, one 
coming from the new cycles of the bubbling geometry and the other given as a surface 
integral on the conformal boundary. It is the former contribution that we manage to 
compute exactly within an ansatz we make. This work does not address the latter 
contribution, which seems to require a holographic renormalization technology beyond 
the one currently available. Indeed, because the new cycles mix non-trivially the AdS^ 
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Figure 1: The Young tableau R is shown rotated and inverted. It consists of g blocks, 
the I-th one of them having nj rows of length Kj. We set = and ^g+i = 

and directions, usual counter-terms in five-dimensional supergravity cannot be used, 
at least in a straightforward way. 

It is however possible to use the identification of the matrix model and supergravity 
data to compare the correlators of the Wilson loop with local operators, namely chiral 
primaries and the energy-momentum tensor. This is reported in a companion paper 

We structure the paper as follows. In Section [2], we study the matrix model for Wil- 
son loops dual to bubbling geometries. We solve the model, obtain its spectral curve, 
and show that it is the hyperelliptic surface that characterizes the bubbling geometry 
dual to the Wilson loop. Section [3] then focuses on the Wilson loop expectation value. 
Using our solution, we compute the Wilson loop expectation value for representations 
that correspond to smooth bubbling geometries. This reproduces the result of [18] in a 
certain limit. We next show that the on-shell supergravity Lagrangian is a total deriva- 
tive and compute the contributions from the new cycles that appear in the bubbling 
geometry. We then conclude the paper by discussing the outlook in Section HI The 
appendices contain details used in the text. 

2 Spectral curves and bubbling geometries 

The expectation value of a circular Wilson loop in A/" = 4 super Yang-Mills is captured 
by a Gaussian matrix model [T2l[T3l[Tl]. This was originally proposed for half-BPS 
loops in the fundamental representation (which are dual to fundamental strings in the 
bulk), but the conjecture has later been extended and checked to hold also for circular 
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loops in arbitrary representations R of the gauge group [T9],[20l[2T],[221[23] and for some 
loops preserving reduced amounts of supersymmetry [2ll[25l[26l[27l[2H]- The precise 
statement is that the Wilson loop expectation value for the U{N) theory is given by 

iWR)u^r,) = [dM] exp (^-^TrM^^ Tr«e^. (7) 

Here M is an hermitian matrix and the partition function Z of the matrix model is 
defined as the integral without the insertion Tr^e^. We use the standard hermitian 
measure [dM]. In the absence of operator insertions, the eigenvalues are distributed in 
the large limit according to the Wigner semi-circle lawjj 

To make better contact with the supergravity solution, it turns out to be more 
convenient to follow the procedure delineated in [11] and decompose M in g + 1 sub- 
blocks M(-f) of size rij X rij. The expectation value of the loop is then given by several 
Gaussian matrix integrals correlated by interactions between the sub-blocks: 

= i / n[«/<"le-¥ E,T.M..,.,.,™"- TT (M<'. « 1 - 1 « A^f 

'J T 1 T ^ J 



7=1 KJ - ^ 



The eigenvalues of M*^^) for fixed I are distributed along some interval [627,627-1]. In 
the following, we drop the exponential interactions by replacing (1 — e~^'^' ® e^'^'^') 
with 1. This is a consistent approximation in the limit 

A»l, gl.^ni = 0{X), g^y^iKj - Kj+,) = 0{\'/^) , (9) 



because 627-1 - 627 = 0{^/9ym^) '^27 - 627+1 = C>{gYM{Ki - Kj+i)) as one can 
see from the saddle point equations below. 

Going to the eigenvalue basis, the matrix model in ([8]) becomes (here i = 1, . . . ,ni 
labels the eigenvalues of the I-th sub-block) 



, (I) (J)] ^ 

I m; — rrij 



(10) 



We have introduced a linear ordering in the set of all the eigenvalues so that the last 
product is taken over distinct pairs of eigenvalues. The saddle point equations are 



4:N (J) 1 



Pedagogical references on general matrix models include 
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By defining the resolvent 



the eqs. ( ITTj) can be written, for x G [627, 627-1], as 

-Ax + glj,iKi + uj+{x)+u.{x) = {), (13) 
where ui±{x) = uj{x ± ie). 

2.1 A hyperelliptic surface as the spectral curve 

By difi^erentiating eq. ( |T3|) . one can see that tu^ = 4 — c<j^, so that the combination 

w'(4 - cu') (14) 

is invariant when crossing the cut. Let us now consider the behavior of this expression 
close to a branch point, say Ci. The eigenvalues are expected to produce square root 
branch cuts. Since uj{z) satisfies eq. f|T3l) . locally it is given by 

uj r^2z- ^gyM^i + cy/z - ei, (15) 
where c is some constant. Then 

u'jA -Lu')r^ (2+ / ] (2 ^ = 4 - , , . (16) 

The same behavior is found for every branch point ef. 

C 

u;'(4 — cij') ~ — as 2; — > Cj , Cj = const., (17) 

so the combination 

2g+2 

z - e,; 



^'(4 

1=1 



is regular everywhere on the complex plane. The first term behaves as 0{l/z'^) for 
large z by the definition of uj. Thus the combination must vanish everywhere and, in 
addition, the second term has to be of the form 



with f2g{z) a polynomial of degree 2g and 



23+2 

H2gUz)^ X{{z-e,). (20) 



i=l 

The solution to the quadratic equation 



u:\A-u') = Jd^ (21) 
-"29+2(2;) 



is then 



a;- = 2-W4--^#-^2-2^ g±i^^ . (22) 
^2,+2(^) 

Here we have selected the negative sign in front of the square root to guarantee the 
correct behavior for z —>■ 00. In introducing the monic polynomial ag^i{z) = z^^^ + . . ., 
we noted that -^23+2 — /2g/4 has to be a perfect square, so that the only singularities 
of uj' are the branch points Cj. 

We can geometrically interpret eq. fl2^ by saying that the resolvent uj{z) is the 
indefinite integral 



uj{z) = a (23) 
J 00 

of a meromorphic 1-form 

a=(2-2^^±M=] dz (24) 
on the hyperelliptic curve defined by 

= H2g+2{z) . (25) 
The only singularity of the 1-form a is the double pole at 2; = 00 on the second sheet. 



2.2 Parameters and constraints 

The parameters in the definition of the spectral curve and the one-form are the 3g + 3 
coefficients of the two monic polynomials a^+i = a and i/25+2 = H. Let us study the 
constraints that determine these parameters. 

The constraints are most concisely expressed in terms of period integrals, so let us 
introduce the A- and -B-cycles of the hyperelliptic surface in the standard way (see fig. 
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Figure 2: The A- and 5-cycles of the hyperelhptic surface S of genus g = 2. 



[2]): the cycle Aj (with / = 1, ...,(? + 1) circles the J-th cut [e2i, €21-1] clockwise. Only 
the first g of the A-cycles are independent, since Ag+i = —Ai — . . . — Ag. The cycle Bj 
(with I = 1, . . . , g) goes through the I-th and the {g + l)-th cuts and has intersection 
numbers #(A/ fl Bj) = 6ij for J = 1,2, . . . , g. 

1. The first g + 1 constraints come from the requirement that the resolvent uj{z) 
should be single-valued on the physical sheet. Since it is obtained by integrating 
the one-form a, we need that 

(f a = 0, I = l,...,g + 1. (26) 

JAi 

These g + 1 constraints are all independent: even though Ai is a trivial 

cycle in homology, the condition J^-, a = applied to (^M is non-trivial and 
ensures that no logarithmic term appears in the expansion of u around z = 00. 

2. According to the saddle point equations ( fT3l) . the value of u along the cycle Bj 
goes from uj to 4z — uj in passing through the {g + l)-th cut from the first to the 
second sheet (recall that Kg^i = 0), and then from Az — uj to uj + gy^Ki when 
coming back to the first sheet across the J-th cut. In terms of the one-form a, 
we get g conditions 

/ a = glMKn 1 = 1,. ..,g. (27) 
Jbj 

3. Since the I-th cut contains nj eigenvalues, the definition (fT2l) implies the following 
g + 1 conditions 

® LU dz = —2nigYj^ni , I = l,...,g + 1. (28) 
Jai 

The integral should be performed on the first sheet. 
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4. The 3(7 + 2 conditions above determine ag+i(z) and H2g+2{z) up to a shift of z. 
The last condition that fixes this ambiguity is 

^{(^23+2) = 2e2g+2 , (29) 
which follows form (fT3|) recalling that Kg^i = 0. 

We check now that uj{z) given by fl23|) together with the constraints fl26|) - fl29|) au- 
tomatically satisfies the saddle point equations (1131) . For this we need to evaluate uj 
just above and below each branch cut [e2/,e2/_i]. Since we know the value of u at 
z = 623+2, we only need to integrate a from 62^+2 to 621 along an arbitrary path on 
the first sheet, and then from 621 to x ±ie with x G [621, €21-1] along the cut. The key 
points are that 

—Lf=dz = gl^,Kj , (30) 

as follows from the condition (1271) . and that 

^/H{x + ie) = -^/H{x-ie) (31) 
on the cut. For x G [621, e2/-i] we have 



uj+{x) + u^{x) = 2w(e2g+2) + 2 / ( 2 - 2—= ) dz 



' 620+2 



a 



+ 1 (2-2^ELUx'+/ 2-2^^\dx' 



e2i,x\+ie 

2 



= 4e2<;+2 + 4(62/ - 623+2) - 9ymKi + ^(^ - ^2/) 
= Ax- gl^Kj , (32) 

so we see that the saddle point equations (fT3!) are satisfied. Thus at this point we have 
found the exact solution of the matrix model (fTO!) in the large N limit. 



2.3 Comparison 

What remains to be shown is that the spectral curve fl25l) is the hyperelliptic surface 
that appears as part of the bubbling solution for a Wilson loop [?]. 

The bubbling geometry is a warped product of AdS2 x S'^ x and a half-plane, 
as we have mentioned in the introduction. This half-plane is taken to be the lower 
half-plane in one sheet of the hyperelliptic surface given by 

2g+l 

s'=l[{u-e,). (33) 

i=l 
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The branch points of the surface are at m = Cj (with i = 1,.. .,2(7 + 1) and u = eo = 
e2g+2 = oo- (Notation changed from jTj: e*'^'''''^ = e^'^'''^.) The constant Uq and the 
branch points Cj are all real and ordered as follows: 

e2g+l < e2g < . . . < ei < Uq. (34) 

Though the RHS of fl33|l is a polynomial of degree 2g + l instead of 2g + 2, the equation 
can be transformed to the form fl2Sl) by a Mobius transformation on u. 

All the supergravity fields are expressed in terms of two holomorphic functions A 
and S on E given by 

3A = ^<">f' (35) 
B = -t^—. (36) 

U — Uq 

The polynomial P{u) has real coefficients and is of degree g + 1. The real part of A 
must vanish on [e2i+i, 621] to ensure regularity of the solution, so there are constraints 

r dA = 0, I = l,...,g + 1. (37) 

The branch cuts [e2i-i, 627-2] represent three-cycles of topology that arise from the 
geometric transition of D5-branes, so they carry RR three-form fluxes. Since each 
column in the Young tableau R represents a D5-brane [201 l3T]. the flux carried by the 
I-th cycle is proportional to Kj — Kj^i, the number of columns in the I-th block: 

8ni f dA+c.c.= f F^f = 47i\Ki - Ki+i)a' (38) 

[627-1,627-2] Js^ 

ioT I = 1, . . . , g. Similarly, the segment [627,627--!] represents a five-cycle of topology 
5*^ that arises from the geometric transition of nj D3-branes [T9l|3T] and carries RR 
five-form fiux. As we show in Appendix [B] 



8ttH I {AdB - BdA) + c.c. = j F(5) = Aii^a'^ni 



(39) 



for J = 1,..., (7 + 1. 

Shifting the imaginary part of A does not affect the physical fields. It is natural to 
fix this ambiguity by requiring that 

lim ^ = 0. (40) 
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The constraints fl371) - P0l) are equivalent to fl26|) - (l29|) . respectively, if we make the 
identification 

8 a' 

uj -2z = i—QsA, B = i—{z - e2g+2)- (41) 

a 4 

Equivalently, we have 

a' r a(z') , , 4 1 

A = i— , ' dz', u-uo = — . 42 

Note that Qym = 4:TTgs. Thus we have showed that the spectral curve of the matrix 
model is precisely the hyperelliptic surface that characterizes the bubbling geometry. 



2.4 SU{N) gauge group 

So far we have focused on the U{N) gauge group case. It is easy to describe what 
changes for a SU{N) theory. First, the Wilson loop expectation value of the gauge 
theory is related to the matrix model by a simple modification of ([7]): 

{WR)suiN) = [dM] exp (^-^TrM^^ Tr^e^^', (43) 

where 

M' = M-l(TrM)l;vxiv (44) 
is the traceless part of M. Since 

Tr^e*^' = e-^T^^Tr«e^, (45) 

the saddle point equation (|T3l) for the J-th cut becomes 

- 4x + g'vMiKj - \R\/N) + co+ix) + co^x) = 0. (46) 

Therefore the resolvents of the U{N) and SU {N) theories are simply related by a shift 
of the argument: 

^su{N){z) = uJu(N){z + |^|/4Ar). (47) 

Equivalently, the eigenvalue distribution is simply shifted by a constant so that the 
average position of the eigenvalues is the origin. The relations between u, z and B 
become 

u; - 2(z + |i?|/4iV) = %-asA, B = t^{z - 62^+2), (48) 

a 4 

where 62^+2 = G2g+2 + |-R|/4A^ is the last branch point in the SU{N) case. 
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3 Wilson loop expectation value 



Given our identification of tlie matrix model and supergravity data, it is natural to 
compare various physical quantities computed on both sides. A companion paper [T7j 
studies the correlators of Wilson loops with local operators, such as chiral primaries 
and the energy-momentum tensor, finding agreement between gauge theory and super- 
gravity analysis. Another natural quantity to compare is the Wilson loop expectation 
value, which we study in this section. On the Yang-Mills side, we compute it using 
our large N solution of the matrix model. We also discuss the supergravity compu- 
tation though we do not complete it in this paperjf] First we prove that the on-shell 
supergravity Lagrangian is always a total derivative. Then we show that the action 
contains contributions from the new cycles of the bubbling geometry and also from the 
boundary of space-time. We compute the first kind of contributions. Issues with the 
second type are discussed in Section HI 



3.1 Wilson loop expectation value from the matrix model 



To the leading order in the saddle point approximation, the normalized Wilson loop 
expectation value is given by 



{Wr) 



-{5niat— 5o) 



(49) 



where iSmat and Sq are the on-shell actions of the Gaussian matrix model with and 
without Wilson loop insertion. We now proceed with computing these actions. 

Again, we begin with the case of a U{N) gauge group. The on-shell value of the 
matrix model action is given by 



-5„ 



E 



2N 



m. 



(I) 



dxp{x 



2N 



{I,i)<{J,j) 



(I) (J) 



[e2l.e2/-i] 



A 



-x^ + Kjx 



where the eigenvalue density 



lA 



dxdy p{x)p{y) log \x - y\, 

(50) 

(51) 



''The computation of the expectation value of a loop dual to D3 and D5 branes [3T1[35] has been 
performed in [l9lf20ll2Tl[22] . both using the matrix model and the DBI action. 
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is related to the resolvent by 



z / pi x] 

p(x) = - — -{ujjf_{x) — u!^{x)) , uj{z) = A / dx- 



27iX 



Z — X 



(52) 



In the limit in which the cuts are well separated, the last term in fISUl) can be dropped, 
and by using the eigenvalue density p(x) given by 



5^ ^5 (x - gl,,K,/A) 



(53) 



we easily reproduce the results of [T8] . 

The expression ( l50|) may be enough for comparison with supergravity although we 
do not see how the double integral can appear in gravity. We now rewrite fISUl) in a 
form that involves no double integral. First, let us use the density and a principal value 
integral to re-express (fTTj) : 



- 4x + qImKi + 2AP / dyp{y)^— = for x G [e2i-i, ea/]. (54) 
This equation can be integrated to yield 

- 2^ + + 2A / log I. - y| ^ 2,?.„o for . e e.]. (55) 



where C/ is an integration constant. The on-shell matrix action is then 



-"mat 



9+1 ^ 

dxp{x) 

1=1 •' [1^21,^21-1] 



^ 2 

~X + -A/X + C/ 
A 2 



(56) 



One expression for the Wilson loop expectation value that does not involve a double 
integral or cj is obtained by using (l55l) with x = 621-1 and x = 62/: 



log{WR)u(N) 

9+1 ^ 

/ dxp{x 



N . 1 
A 



K A 621-1 + 62/) 



/=1 [e27,e2/-l] 

+ ^ X] '^J log(e2J_l -X){X- 62J 



log VA + 3/4 + log 2, 



(57) 



where we used 

5o = (-logyA + 3/4 + log2) (5J 
that follows from the density po{x) = (I/ttA) VA — x^ for Wigner's distribution. 
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For the SU{N) theory, we simply replace Kj by Kj — \R\/N: 





-logVA + 3/4 + log2. 



(59) 



In this formula p{x) and Cj are the density and the branch points in the SU{N) case, 
and we have used the fact that the average eigenvalue vanishes to remove a shift of Kj 
in the second term inside the bracket. 

3.2 Wilson loop expectation value from supergravity 

Let us now turn to supergravity. The solution in [7J is for an infinite straight line along 
the Lorentzian time, whereas the matrix model model computation is appropriate for 
a circle in Euclidean signature. This is not a problem, since both the straight line 
and the circle preserve the same isometry S0{2, 1) x S0{3) x 5*0(5) (albeit reahzed 
differently in the two cases). We can then extend the solution of to the circular case 
via a Wick rotation, considering a fibration with the Euclidean factor HI2, rather than 
AdS2- This difference will not play any significant role in our analysis, so that we shall 
consider for simplicity the Lorentzian signature. The Wilson loop expectation value is 
then given by {Wr) = exp(— 5^;) after the Wick rotation that identifies —Se with iSi, 
where Se and Sl are the Euclidean and Lorentzian on-shell actions. 

3.2.1 The on-shell Lagrangian is a total derivative 

We begin our discussion of the supergravity action by showing that the on-shell La- 
grangian density always has to be a total derivative, if it is a homogeneous function 
of the fields of non-zero degree. It seems well-known that the supergravity Lagrangian 
is a total derivative if the equations of motion are satisfied, though we do not know a 
reference that makes the general statement explicitly. 

The argument is simple. Suppose the Lagrangian C{4>) depends on the fields 0* 
and their derivatives. There can be second or higher derivatives. When we take the 
variation of C with respect to arbitrary changes 5(j)'^, in general we get terms that contain 
derivatives of 6(f)\ By definition, the equations of motion Si{(j)) = are obtained by 
rewriting 6C as 




(60) 
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where T>i is the total derivative term that is hnear in S(j)^. If the Lagrangian is homo- 
geneous, there are (usually integers) numbers ric and Ui such that 

£(fi"»0O = Q"^C{(f)') (61) 

for any constant Q. We call Ui the dimensions of the fields. By choosing = 1 + e so 
that 50* = eni(f)^, we find that 

enz:/:(0) = J]eni0*£:i(0) + ^A(en,0^0). (62) 

i i 

If the equations of motion are satisfied, the Lagrangian is a total derivative: 



We now apply the above consideration to the type IIB superg ravity actioifl 
J ^ ^\ 2 (Imr)' ; 

+ J (^-^M,,H^,^ A ^Hl,^ - 4F(5) A *F(5) - e,,C(4) A H^,^ A /ffg) j . (64) 

The action is written essentially in the convention of [M] and contains various combi- 
nations of the fields: 

T = C(o) + ze-'^ , (Mab) = diag(e-^ e^) , 

F(5) = ^^(4) + A dBl,^, (65) 

where H^^^ = dB'^2) a = NS, RR. First note that the action is homogeneous of 
degree 8 if we assign dimension 2 to the metric qmn and p to all p-form fields (scalars 
are zero-forms). So our argument applies. Since the scalars have vanishing dimensions, 
we can ignore their variations. Then under arbitrary variations of the fields except the 
scalars, the action changes as 

2n^5S = J d^\^V^\v''5gMN-g'''^VM5gpQ) 

+ jd(^- MabSBl^) A ^Hl^ - 2e,,C(4) A SBl^^ A i/fg) 

- 85C(4) A ^F(5) + eabSBl^) A B\^^ A *F(5)^ (66) 



®Sclf-duality of the five-form, F^^-^ — *F(5), does not follow from this action, but has to be imposed 
by hand. One can consider other actions where self-duality is implemented with an auxiliary field. In 
the case [33] we looked at, the on-shell value does not seem to change. 
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up to terms that vanish on-shell. By setting 

5gMN = 2egMN, 5B^^^ = 2eB^^^, 5C(4) = 4eC(4) (67) 
and using SS = 8eS, we conclude that the on-shell action is given by 

2k^S = j d (^-lM,,5f2) A ^iffg) - ^e„,q4) A Bl^) A iffg) - 4C(4) A *F(5)^ . (68) 

We thus see that we only need the two- and four-form fields to compute this part of 
the action. Note that so far we have not committed to any particular solution. 



3.2.2 Contributions from new cycles 

In the solution of |7], the NS two- form is along the AdS2 directions while the RR 
two-form along the 5*^ directions. The RR four-form has two components, one in the 
AdS2 X S'^ and the other in the S"^ directions. One has then 

dNS _ L -01 rRR _ l s23 n — ^ 20123 , ■ -4567 ( (!(\\ 

where e°^, e^^, and e^^^^ are the volume forms of 74^5*2, 5^, and S"^, respectively, all 
normalized to unit radius. Note that 61, 62, Ji, and j2 are real functions on S. Recall 
now that the S'^ and S"^ radii vanish on segments of the real axis of E. Thus e^^ and 
^4567 g^j^g j^Q^ globally defined forms in the ten-dimensional space-time, while e^^ is. 
This implies that the Chern-Simons term in flM|) is not globally defined. We can make 
it globally defined by adding further total derivative terms 

2^S, = J d (2^(4) A 5(f A i/gf - ^Bf,^ A 5(f A d {B^f A 5gf )) (70) 

so that the new Chern-Simons term in 2/t^(5 -|- Si) = is 

J 2F(5) A A H^f. (71) 

The on-shell action is then given by 

2K'5b,ik = Jd (^-^Ma,B^,^ A ^iffg) - A B^f A dC^,)^ , (72) 

where we took into account flB^ . 

Since some of the forms in fl72l) are not globally defined, we need caution in applying 
the Poincare lemma. Some terms in fl72l) are contributions of the non-trivial cycles in 
the bubbling geometry, while the rest are from the boundary of space-time. We focus 
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on the former contributions. The latter should be combined with counter-terms we do 
not discuss in the present work. 

With our ansatz, the Hodge duals of the three-forms are given by 



^Fgf = Ae~, (73) 

Ji 



NS _ J2J4 , TT . -234567 



^F^f = h^^dhAe'''^'''. (74) 

J2 



5) ~ 

We have by the Poincare lemma 

2K25Mik = ^ (-i^"""^^! * - 'a''^^'^^^ * "^^2 ~ hhdj^ . (75) 

By 9S we denote the real axis as well as a large semi-circle on the lower half-plane. We 
cannot meaningfully separate contributions from the two components of dT, because 
adding an exact form in the integrand of fl75p mixes them. In 0751) we have made the 
important assumption that the volume of AdS2 is regularized in a way independent of 
the position on 9S. We have denoted the volume of AdS2 x 5*^ x 5*^ by V. In a more 
complete calculation of the on-shell action, this assumption may need to be modified. 

In Appendix [21 we study how various quantities in fl75p behave in the asymptotic 
region z —>■ oo. If we choose the coordinate to be the spectral parameter z in the SU {N) 
case, both bi and 62 vanish as 2; — >■ cxo while j2 remains finite. Thus the contribution 
from the semi-circle in this parametrization vanishes. 

On the real axis, the first term in fl75|) never contributes because it contains positive 
powers of radii of the two spheres and always vanishes. The remaining two terms nicely 
combine to give 

00 / 1 ^2/4 

4 



V ^2 76'^^ ^ db2 + hd32 ■ (76) 



The sign change from f[75l) is due to the natural direction for integration. We observe 
that /4 vanishes on regions of the real axis where shrinks to zero size. In fact, j2 is 
constant there because otherwise F5 that contains dj2 A e'^^^'^ would be ill-defined. On 
the other hand, 62 is constant on regions where S'^ shrinks for a similar reason and, as 
we recall in Appendix lAl 62 = — 4Im^. Since ^(62^+2) = by (140!) . the flux condition 
(!38|) determines these constants to be 

b2 = 2n:a'Kj on [621,621-1]. (77) 

Thus we can write fl7^ as 

<?+l /. / .2/4 



V / * db2 + Abidj2 

1=1 ^ J[e2l,e2i-i] \ J 2 



(7J 
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equation of motion for can be written as 



The physical meaning of the integrand in flTHl) can be understood as follows. The 

^2) 

where 



dHj^^ = 0, (79) 



H^f ^ ^ H^f + 4i?(f A - A 5gf A H^f. (80) 

It is easy to see that the integrand in (!78l) is proportional to the component of 
along the AdS2 x S'^ direction. The seven-form is to be regarded as the field strength 
of the six-form potential = dC^. By the symmetry of 74^5*2 x S*^ x 5*^, we can 
write 

where e°^'^^^'' is the volume form of unit AdS2 x S^. Then by definition 



£2 £4 

^ill^ * dh2 + 46idj2 = rf&4. (82) 

J2 



e 



Thus the integrand in (ITS]) is ^64. 



One can express the LHS of (1821) in terms of A and z using the known expressions 
for fields summarized in Appendix [XI It is in fact possible to integrate the equation: 

J_ _ 2{z - z){A + Af - {z^ - z^){A + A){d,A + M) 
a'2 ^ ~ 2{d,A - d-,A) 

^-h^z^ -z^){A- A)-^ J dzzA-6 j dzzA, (83) 

where the last two terms involve indefinite integrals. One can check that (!82|) is satisfied 
by this solution. On the real axis where z = z, b4 reduces to 



64 = -6a'2 j dzzA + c.c. . (84) 

Thus 

&4(e2/-i) - &4(e2/) = -evr^a'^A^ / dxp{x)x . (85) 

[e2l,e2l-l] 

By collecting everything together, (175!) becomes 

1=1 he 



2K^Sbnik/V = --Tt^a'^N} ^ Ki I dxp{x)x. (86) 



9 / ■> I 

1=1 -'[e2l,e2l-i] 

This is the contribution from the bulk, in particular from the cycles that have grown 
in the bubbling geometry. This is not the complete story, since the volume V should 
be regularized and counter-terms on the boundary should be added. We see indeed 
that fISB]) seems to account only for special terms in the matrix model action in flS^ . 



18 



4 Conclusion 



The main achievement of this paper is the large solution of the matrix model that 
governs circular BPS Wilson loops at strong coupling. We determined the eigenvalue 
distribution for an arbitrary representation in terms of geometric data on the spectral 
curve. The spectral curve was then identified with the hyperelliptic surface S that was 
found in [7] to characterize the bubbling geometry for the Wilson loop. 

The identification of the hyperelliptic surface S as a spectral curve is important 
for two reasons. First, one can view this as an example of emergent geometry. The 
matrix model is a reduction of the four- dimensional gauge theory [H] and the geometry 
emerges out of the dynamics of the eigenvalues. 

Second, the identification provides the precise dictionary between field theory and 
gravity. Indeed it serves as the basis for the matching of physical quantities computed 
on both sides. A successful example of matching is reported in [T7|, where the corre- 
lators of the Wilson loop with chiral primaries and the energy-momentum tensor are 
computed. 

It should also be possible to match the computations of the Wilson loop expectation 
value. Given our solution of the matrix model, we were able to compute the Wilson loop 
expectation value quite easily. On the other hand, the computation of the expectation 
value in supergravity is unfinished. Such computation should involve two non-trivial 
tasks. One is to properly take into account the new cycles that appear in the bubbling 
geometry. In the present work, we developed techniques to perform this task. The other 
task is to regulate the infinite volume of the ten-dimensional space-time and to add 
proper counter-terms. Usual five- dimensional counter-terms do not suffice, because the 
bubbling geometry mixes the AdS^ and directions in a topologically non-trivial way0 
Construction of the counter-terms is a worthwhile open problem that has applications 
to other observables such as surface operators [37]. 

Acknowledgments 

We are happy to thank David Berenstein, Eric D'Hoker, Sean Hartnoll, Gary Horowitz, 
Don Marolf, David Mateos, Joe Polchinski, and especially Kostas Skenderis for fruitful 
discussions. We are also grateful to Jaume Gomis and Shunji Matsuura for collab- 
oration on a related project. This research is supported by the NSF grants PIIY05- 
51164 and PIIY04-56556, and by the Department of Energy under Contract DE-FG02- 

similar problem, related to the difficulties of formulating higher-dimensional counter-terms [35], 
was already encountered by one of the present authors in the context of bubbling geometries [36] . 
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A Details on the bubbling geometry 



The solution to the BPS equations can be expressed in terms of two holomorphic 
functions A and B on the lower half-plane. Let us define four harmonic functions 
hi, hi, /i2, and /12 by 



A=-{hi-thi), B 



^{h2-ih2) 



^7) 



In fact, all the physical fields except the form fields can be written in terms of hi and 
/i2 alone. The field strengths of the form fields are also given in terms of hi and h2. 
The dual harmonic functions hi and h2 only appear in the potentials [7j. 

It is useful to define the following shorthand notations 



V = d^hidtoh2 - du,hid^h2 

|2 7,2t 



W = d^hidu,h2 + du,hid^h2 

2 z,2t 



Ni = 2hih2\d^hi\' - hfW, N2 = 2hih2\d^h2 
where w is an arbitrary complex coordinate on S. Then we have 



(88) 



12 

iVi' 



p 



w 



h\h\ 



ft = -4e^/^l^ , = 4e-^/^^^ , f! 

while the relevant components of the two- and four-form fields ( l69l) are 

hih^V 



No 



as given in [7], and 



bi = -2^%^ _ 2h2 , b2 = -2i-^^ + 2hi 

iVi 1^2 

V - 3 - 

32 = ^^1^2^ + - - 2(^1^2 - hih2) 



^9) 



(90) 



(91) 



as we show in Appendix [Bl The holomorphic function C is defined implicitly by 

a„C = Ad,„B - Bd,„A . 



(92) 



The behavior of various quantities near the real axis {y = 0) was studied in [7j: 



Intervals 


Vanishing fiber 


hi 


dyhi 


h2 


V 


W Ni 


N2 


[e2/,e2/-i] 




0{1) 


0{y) 


0{y) 


0{l) 


0{y) Oiy) 


0{y) 


others 




0{y) 


0{l) 


0{y) 


0{y) 


0{1) 0{y') 


0{y') 



It follows that 62 = 2/^1 = — 4Im^ on [e2i, 627-1] 
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B An explicit expression for the four-form 



The component j2 of the RR four- form C(4) (l69l) is not given exphcitly in [7], but can 
be obtained along the hues of the similar computation in Section 9.9 of [38]. We use 
the notation of these papers. 

The derivative of j2 admits an expressioii^ 

dj2 = -ift {pfwdw - pfujdw) , (93) 

where from eqs. (5.24) and (6.1) of [7J and from the relation pp^ = dw4> (0 = (p/2) one 
has 

2pf^ = log ^+(^-^] 9^0. (94) 

Using that 



a \aa (3 (3 J 



a = ^ '^y cosh{(j) + X) , /? = ^^/sinh(0 + A) , (95) 



it becomes 



r, r dnj(p + dyjX 2cosh(A-A) 
smh(20 + 2A) | smh(20 + 2A) | 

The warp factor is given by /4 = z/(q;/3 + /3a) (with u = ±1), so that 

2 — 2 

/4 = 1^ (sinh(20 + A + A) - I sinh(20 + 2A)|)^ . (97) 

One can now change the variables from (real) and A (holomorphic) to the real vari- 
ables /i and -(9 defined by 

, T . sinh(20 + 2A) 

A-A = «/i, e = ■ T,/o , , ' 

smh(20 + 2A) 

from which also follows 

|sinh(20 + 2A)P = -— ^-..^ |sinh(20 + 2A)| 

' ^' 4 sin(7? + /i) sin(t9 - /i) ' sinh(20 + 2A) ^ ^ 



and 



sin2/i(9^„^? z sin2i9 9^/x 

4 sm^t/ + yUj sm[v — fi) 2 4 sm^t/ + /ij sm^t/ — yUj 



Using complex coordinates on E, the frames become — e^" + e™, = —1(6™ — e™), with 
: pdw and e™ = pdw. 
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Using eq. (7.4) of [7] one has 
1 



2/)'^ cos^ fi - 



e * (— sin2yU(9^'(9 — ze^* + z cos2yU(9^/i) — 
-2 cos /i (— sin 2/i S^-i? + ie'"^^"^ dujjJi — i cos 2/i . (101) 



In terms of = ^'^Z^, the expression above becomes 



2i 



(sin2/i)2 L 



— sin 2/i (^'ip dw-ip — ip^ dyjip /ip) — ip'^ + cos 2/i'?/'^ dwP+ 
+2 cos /X sin 2/i (-^ S^u?/^ — ip d^jip) — 2 cos /i <9^/i / "^H- 

(102) 



+2 cos /i cos 2/i c^iu/i 
and finally, using the equation of motion 



dw^p = cot /i V' H f^to/^ 

sm /i 

to eliminate the pieces with more than 2 ip and/or ip, 



(103) 



sin 2/i 

This can be almost written total derivative 



^pd^ ^ ^jj^^jj^^ ^ 2 cos 2/i ^2 g ^ 

sin 2/i (sin2/i)2 (sin2/i)2 

2 cos /i , - ^ , , „ 2 cos /i cos 2/i - ^ 

[i) d^ip - i) d^ip) H . ■ _ — V'V' OuiP 



sin2/i)2 

ip'^d^H 



2p/u,/4 = ^^^^ - 3z .2 

sm/i sm 2/i / sm /i 



(104) 



(105) 



using again the equation of motion for ip. Using the equation of motion for ip^ eq. (7.7) 
of [7], the expression for k in eqs. (7.8) and (7.13) and the last equation in (7.14), the 
last term in the formula above becomes 

ip'^d^H 



9 — dw\ip cot p + ihiC 
sin /i 



+ 2t{hM2 - h2du,hi) . (106) 



Then 



and one has 



pfwft = -idy,{hih2 tan fi) - 3(/ii9^/i2 - h2dyjhi) 



- 3 ~ 

j2 = -/ii/i2tan/i + 32(C - C) - -(/ii/i2 - ^^1/^2) • 



(107) 



(lOJ 



Using ( l92l) together with the relations /i = — ^(A — A) and e^^ = dy,hi/dwh2, one can 
rewrite this as (1911) . 
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C Asymptotic behavior 



Let us study the asymptotic forms of physical fields in the region z ^ oo. We use the 
SU{N) identification (HHj) of the matrix model and geometry data. 

From the definition f[T^ . u behaves in the asymptotic region of S as 

u;iz) = - + o(\]. (109) 



z \ z^ 



The order 0{z ^) term vanishes in the SU{N) case. Using the formulas in Appendix 
[K\ we find the asymptotic forms of various fields: 

e^ = e2^ = ^?, + 0(r-4), (110) 
r + 0(l/r), f^=[jx) ^ + (111) 



/4=(^^J \sme\ + 0{l/r'), (112) 

6i = 0(l/r), 62 = 0(l/r), (113) 
^ _^-A[12.-8sm(2.)+sin(4.)] _^ ^^^^^.^^ ^^^^^ 

64 = C(l/r). (115) 

Here we introduced polar coordinates z = re*^ with —it < 9 < 0. Note that the 
metric ([1]) is written in the Einstein frame where the AdS radius is {a''^X/ gsY^^ in our 
convention (11 101) for the dilaton. The subleading terms depend on the representation 
R and can be easily calculated in terms of the moments of the eigenvalue distribution. 



References 

[1] D. Berenstein, "Large N BPS states and emergent quantum gravity," JHEP 01 
(2006) 125, [arXTv • hep-th/0507203[ 

[2] D. Berenstein, "A strong coupling expansion for N=4 SYM theory and other 
SOFT'S," [arXiv: 0804 .0383 [hep-th]] 



[3] D. Berenstein, "A toy model for the AdS/CFT correspondence," JHEP 07 



(2004) 018, arXiv : hep-th/0403110, 



[4] H. Lin, O. Lunin, and J. M. Maldacena, "Bubbling AdS space and 1/2 BPS 



geometries," JHEP 10 (2004) 025, arXiv:hep-th/0409174 



23 



S. Yamaguchi, "Bubbling geometries for half BPS Wilson lines," 



Int. J. Mod. Phys. A22 (2007) 1353-1374, larTiv : hep-th/0601089 



O. Lunin, "On gravitational description of Wilson lines," JHEP 06 (2006) 026, 
larXiv : hep-th/0604133| 

E. D'Hoker, J. Estes, and M. Gutperle, "Gravity duals of half-BPS Wilson 
loops," JHEP 06 (2007) 063, . arXiv: 0705. 1004 [hep-th] , 

J. Gomis and S. Matsuura, "Bubbling surface operators and S-duality," JHEP 



06 (2007) 025, arXiv: 0704 .1657 [hep-th] 



J. Gomis and T. Okuda, "D-branes as a Bubbling Calabi-Yau," 



JHEP 07 (2007) 005 | [ arXi v : 0704 . 3080 [hep-th] [ 

J. Gomis and T. Okuda, "Wilson loops, geometric transitions and bubbling 
Calabi- Yau's," JHEP 02 (2007) 083, arXiv:hep-th/0612190, 

N. Halmagyi and T. Okuda, "Bubbling Calabi-Yau geometry from matrix 



models," \JHEP 03 (2008) 028| |arXiv: 0711. 1870 [hep-th] 



J. K. Erickson, G. W. Semenoff, and K. Zarembo, "Wilson loops in N = 4 



supersymmetric Yang-Mills theory," NucL Phys. B582 (2000) 155-175 



arXiv : hep-th/0603055[ 



N. Drukker and D. J. Gross, "An exact prediction of N = 4 SUSYM theory for 



string theory," |J. Math. Phys. 42 (2001) 2896-2914, arXiv:hep-th/0010274 



V. Pestun, "Localization of gauge theory on a four-sphere and supersymmetric 
Wilson loops," arXiv: 0712. 2824 [hep-th] f 

G. Bonelli and H. Safaai, "On gauge/string correspondence and mirror 
symmetry," |arXiv r0804 . 2^2¥ [hep-th]^ 

N. Berkovits and C. Vafa, "Towards a Worldsheet Derivation of the Maldacena 



Conjecture," JHEP 03 (2008) 031, |arXiv: 0711. 1799 [hep-th] 



J. Gomis, S. Matsuura, T. Okuda, and D. Trancanelli, "Wilson loop correlators 
at strong coupling: from matrices to bubbling geometries," 
[arXiv : 0807 . 3330 [hep-th] , 



J 



T. Okuda, "A prediction for bubbling geometries," \JHEP 01 (2008) 003', 



[arXiv: 0708. 3393 [hep-th] 



24 



[19] N. Drukker and B. Fiol, "All-genus calculation of Wilson loops using D-branes," 



\ JHEP 02 (200 5) 010| |arXiv:hep- th/050 1109| 



[20] S. Yamaguchi, "Wilson loops of anti-symmetric representation and D5- branes, 
JHEP 05 (2006) 037, larXiv : hep-th/0603208 



[21] K. Okuyama and G. W. Semenoff, "Wilson loops in N = 4 SYM and fermion 



droplets," JHEP 06 (2006) 057, ,arXiv:hep-th/0604209 



[22] S. A. HartnoU and S. P. Kumar, "Higher rank Wilson loops from a matrix 
model," JHEP 08 (2006) 026, |arXiv : hep-th/0605027[ 

[23] S. Giombi, R. Ricci, and D. Trancanelli, "Operator product expansion of higher 
rank Wilson loops from D-branes and matrix models," JHEP 10 (2006) 045, 



arXiv : hep-th/060F077 



[24] N. Drukker, "1/4 BPS circular loops, unstable world-sheet instantons and the 
matrix model," JHEP 09 (2006) 004, arXiv : hep-th/0605151 . 

[25] N. Drukker, S. Giombi, R. Ricci, and D. Trancanelli, "On the D3-brane 
description of some 1/4 BPS Wilson loops," JHEP 04 (2007) 008, 
[arXiv:hep-th/0612168f 

[26] N. Drukker, S. Giombi, R. Ricci, and D. Trancanelli, "More supersymmetric 
Wilson loops," ]Phys. R ev. D76 (2007) 107703, a rXiv: 0704. 2237 [ hep-"th]] 

[27] N. Drukker, S. Giombi, R. Ricci, and D. Trancanelli, "Wilson loops: From 



four-dimensional SYM to two-dimensional YM," \Phys. Rev. D77 (2008) 047901 



arXiv: 0707. 2699 [hep-th] 



[28] N. Drukker, S. Giombi, R. Ricci, and D. Trancanelli, "Supersymmetric Wilson 
loops on 53," JHEP 05 (2008) 017, ,arXiv: 0711 .3226 [hep-th], 

[29] P. H. Ginsparg and G. W. Moore, "Lectures on 2-d gravity and 2-d string 



theory," hep-th/9304011 



[30] P. Di Francesco, P. H. Ginsparg, and J. Zinn- Justin, "2-d gravity and random 



matrices," Phys. Rept. 254 (1995) 1-133, |hep-th/9306153 



[31] J. Gomis and F. Passerini, "Holographic Wilson loops," JHEP 08 (2006) 074, 
' arXiv : hep-th/0604007, 

[32] J. Gomis and F. Passerini, "Wilson loops as D3-branes," JHEP 01 (2007) 097, 
|arXiv : hep-th/0612022f 



25 



[33] G. Dall'Agata, K. Lechner, and D. P. Sorokin, "Covariant actions for the bosonic 
sector of D = 10 IIB supergravity," Class. Quant. Grav. 14 (1997) L195-L198, 



arXiv : h ep-th /9707044[ 



[34] E. D'Hoker, J. Estes, and M. Gutperle, "Ten-dimensional supersymmetric Janus 
solutions," Nud. Phys. B757 ( 2006) 79-116| [iLrXiv:hep-th/0603012 



[35] M. Taylor, "Higher-dimensional formulation of counterterms," 
arXiv:hep-th/0110142, 

[36] S. Giombi, M. Kulaxizi, R. Ricci, and D. Trancanelli, "Half-BPS geometries and 
thermodynamics of free fermions," JHEP 01 (2007) 067, 
[arXiv:hep-th/0512101( 

[37] N. Drukker, J. Gomis, and S. Matsuura, "Probing N=4 SYM With Surface 
Operators," |arXiv: 0805 .4199 [hep-th] ] 

[38] E. D'Hoker, J. Estes, and M. Gutperle, "Exact half-BPS Type IIB interface 
solutions I: Local solution and supersymmetric Janus," JHEP 06 (2007) 021, 



arXiv : 0705 . 0022 [hep-th] 



26 



